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Abstract. It is becoming increasingly clear that galaxy bias is compli- 
cated, with evidence indicating that it can be time-dependent, morphology- 
' dependent, scale-dependent, non-linear and non-deterministic. We discuss 

ly-} ■ strategies for overcoming these complications and performing precision cos- 

mological tests with upcoming redshift surveys. 
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1 The Bad News 



The relative distribution of galaxies and mass is of increasing concern 
in cosmology since constraints on cosmological parameters from galaxy 
surveys |L5], [j], ||, § are only as accurate as our understanding of bias. 
Q> . Unfortunately, bias is complicated. The commonly used assumption 

that the matter density fluctuations 5(r) and the galaxy number density 
Qh' fluctuations g(r) obey 

6 : g(r)=bS(r) (1) 

h ; 

"£5 ! f° r some constant b (the bias factor) appears to be violated in a number 

of ways. It has been known since the outset || that b must depend on 
galaxy type. However, there is also evidence that it depends on scale 
(see e.g. [11, |jj and references therein) and on time. Time-evolution 



$— 1 ' of bias is both expected theoretically ||, [l6| and observed Q, with b 

generally dropping towards unity from high values early on. This is due 
to two separate effects: (1) Early galaxies are more biased at birth since 
they tend to be associated with rarer peaks. (2) Old galaxies gradually 
debias as they are gravitationally pulled into the same overdensities 
as the dark matter (see Figure 1). Moreover, on small scales where 
non- linear dynamics is important, one expects nonlinear corrections to 
equation (|l|). Finally, there are good reasons to believe that there is 
no deterministic relation that can replace equation (|l|), but that bias is 
inherently somewhat stochastic [0, 12, M. 



2 The Good News 



Does this mean that galaxy redshift surveys will contribute nothing to 
the quest for accurate measurements of cosmological parameters? Not 
necessarily. Below we will describe an approach that may be able to take 
us beyond the murky marshes of bias and let us do precision cosmology 
without assuming any bias relation in place of equation ([l]) . 

Dekel and Lahav have proposed a framework termed stochastic bi- 
asing which drops the assumption that the galaxy fluctuations g(r) 
is uniquely determined by the matter fluctuations S(r). Instead, g is 
modeled as a function of S plus a random term. Restricting attention 
to second moments, all the information about stochasticity is contained 
in a single new function r(k) (l^, |l6|| . Grouping the fluctuations into a 
two-dimensional vector 



and assuming nothing except translational invariance, its Fourier trans- 
form x(k) = J e _jk ' r x(r)(i 3 r obeys 

(x(k)x(k0t) = (2^(k-k')( £W £M) (3) 

for some 2x2 power spectrum matrix that we will denote P(k). Here 
P is the conventional power spectrum of the mass distribution, P g is 
the power spectrum of the galaxies, and P x is the cross spectrum. It is 
convenient to rewrite this covariance matrix as 

p « = ^ ( KkMk) li m ' ) (4) 

where b = (Pg/P) 1 ^ 2 is the bias factor (the ratio of galaxy and total 
fluctuations) and the new function r = P x / {PP g ) 1 / 2 is the dimension- 
less correlation coefficient between galaxies and matter. Note that both 
b and r generally depend on scale k. The Schwarz inequality shows that 
the special case r = 1 implies the simple deterministic equation ([[]), and 
the converse is of course true as well. 

Since the function r(k) is a cosmologically important quantity, it has 
received much recent attention. Pen has shown how it can be measured 



using redshift space distortions and nonlinear effects [12|, and it has 



recently been computed for a number of theoretical models [13]. Its 



time-evolution has been derived in the linear regime [|16]] and generalized 



to the perturbatively non-linear case [111]. For a fixed galaxy population 
in an in an Einstein-de Sitter universe (O m = 1 and S7a = 0, as in the 
standard cold dark matter model), the relation between the values (6, r) 



at redshift z and the present values (&0) r o) reduces to [16| 



[z 2 -2z(l+z)b r + (l + z) 2 b 2 } 1 / 2 , 
[(1 + z)6 r - z]/b, 



(5) 
(6) 



plotted in Figure 1. Both b and r have now been measured using hy- 
drodynamic simulations M, and other groups are performing related 
studies. 
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Figure 1. The 

evolution of bias (A) 
and correlation (B) is 
shown for 15 models 
that at Einstein-de 
Sitter redshift z = 5 
have bias b = 0.5, 1, and 
2 and correlation 
coefficients r = 0, 0.25, 
0.5, 0.75, and 1. In the 
top panel, r increases 
upward in each 
quintuplet of lines. In 
the bottom panel, bo 
increases downward in 
each triplet of lines. 
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Why is this useful? How can one learn anything about the matter 
density without making assumptions about how it relates to the galaxy 
density? By assuming that the velocities are the same for galaxies 
as for other matter on large scales, i.e., that they are both caused 
by gravity. This means that the characteristic clustering anisotropies 
known as redshift-space distortions ^] will appear even if galaxies 
are completely uncorrelated. Within the stochastic biasing framework, 
Kaiser's classic result (which assumes r = 1) gets generalized to [jl2], |2| 



PobsQt) = Kk)* + 2n°Jb{k)r(k)x + o 



'm x 



PQ*), 



(7) 



where x = (f -k) 2 andf is the line of sight direction (for the special case 
of a volume- limited distant survey). The right-hand side is a quadratic 



polynomial in i, so with a large data set such as the 2-degree Field 
Survey (2dF) or the Sloan Digital Sky Survey (SDSS), one can accu- 
rately measure the three coefficients. For each scale k, equation ([?]) 
thus gives us three independent equations. If the matter density f2 m 
can be accurately measured by other means (using say upcoming CMB 
experiments, supernova la surveys and cluster abundance evolution), 
then these three equations are readily solved for the three unknowns 
P(k), b(k) and r(k). The last two would then provide valuable insights 
about galalaxy formation, whereas P (k) would give us the long sought 
for power spectrum of the underlying matter distribution which could, 
among other things, provide precision constraints on neutrino masses 
H and the Hubble constant j|] . 
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